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Abstract 


One of the major problems in applied mathematics and engineering 
sciences is solving nonlinear equations. In this paper, a family of eight- 
order interval methods for computing rigorous bounds on the simple zeros 
of nonlinear equations is presented. We present the convergence and er- 
ror analysis of the introduced methods. Also, the introduced methods 
are compared with the well-known interval Newton method and interval 
Ostrowski-type methods. Finally, we propose a technique based on the 
combination of the newly introduced approach with the extended interval 
arithmetic to find all of the roots of a nonlinear equation that are located 
in an initial interval. 
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1 Introduction 


The main motivation for this study is to enclose the simple root «* of the 
nonlinear equation 
f(x) (1) 


—d 0, 
by a bounded interval, where f : D C R — R is a real-valued nonlinear 
function on the open interval D. 
Nonlinear problems are of interest to engineers, physicists, and many 
other scientists because most systems are inherently nonlinear in nature. 
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Up to now, many modified methods for solving nonlinear equations have 
been developed to improve the local order of convergence of some classical 
methods, such as Newton, Chebyshev, Potra-Ptak, and Ostrowski methods; 
see [19, 18, 7, 8, 6, 3, 4, 9, 10, 14, 2, 23, 5, 13). 


An optimal eight-order method for solving nonlinear equation (1) pro- 
posed by Bi, Ren, and Wu [2] that is based on King family [14], is given 
by 

Yn = Ty — fee , 
= 2f(@n)—-f(yn) f(yn) 
2m = Un — DFE =BF Ya) Pen) (2) 


In+1 = %n — (Mn) GX 


where [ln = oe and his a real-valued function with h(0) = 1, h’(0) = 2 and 
|h’"(0)| < co. Iterative method (2) with eight-order of convergence is very fast 
compared with many other methods. Solving the problems in floating-point 
arithmetic is inevitably associated with round-off errors, and so the obtained 
solution to the problem is accompanied by some errors. Interval analysis is a 
tool for bounding the errors and providing rigorous bounds on the solution to 
the problems. The interval extension of the Newton method with quadratic 
convergence [24, 16], the interval extensions of the Ostrowski method and 
modified Ostrowski method, respectively, with fourth-order and sixth-order 
of convergence [1, 11], and the interval extension of the n-step Traub method 
with (nm + 1)-order of convergence [21], are examples of the interval methods 
that give rigorous bounds on the solution to the nonlinear equations. 


In this work, we present an interval extension of (2), which has an eight- 
order of convergence and gives rigorous and outstanding results, that is, in- 
terval enclosures with sharp bounds that contain the exact solution. Also, we 
introduce a technique based on combining the new method and the extended 
interval arithmetic for enclosing all simple roots that are located in an initial 
interval. In contrast, many root-finding methods can only find one root of 
the function in the given initial interval. 


Here, we use boldface letters to denote intervals. The set of real intervals 
is denoted by IR = {x = [x,x] : x < x}. The midpoint and width of an 
interval number x = [x,x] are defined by m(x) = as and w(x) = X—x, 
respectively. The absolute value of x is |x| = max{|z| : « © x}. The interval 
extension of real-valued function g is denoted by its corresponding uppercase 


and bold letter G. 


2 Description of the methods 


Many modified methods for solving nonlinear equation (1) with a high-order 
of convergence are based on the well-known Newton method. So, we first 
give a brief description of the interval Newton method. 
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2.1 Interval Newton method 


The idea of the interval Newton method for the first time was discussed 
in [24, 16]. Suppose that the real differentiable function f in (1) has the 
inclusion of monotonic interval extension F’(x) of its derivative f’(a) and 
that xo is an initial point. Then the interval Newton method is 


Flom) yg 


P(x) ma 0; 1 ass (3) 


Xn+1 = {m(x») _ 
Recursive relation (3) produces a sequence {x,,} of interval numbers. If the 
initial interval x9 contains a zero «* of f(a) and 0 ¢ F’(xo), then all iterates 
contain «* and the method converges to x”*. 


Theorem 1. [17] Let f be a real rational function of a single real variable 
x with rational extensions F and F’ of f and f’, respectively, such that f 
has a simple zero y in an interval [21,22] for which F'([@1,x2]) is defined and 
F’((x1,%2]) is defined and does not contain zero. Then there is an interval 
Xo C [21, £2] containing y and a positive real number K such that 


W(Xn41) < K (w(Xa))’, 


therein {x,,} is the produced interval sequence by (3). 


2.2 Main results and convergence analysis 


In this subsection, a new interval method is introduced to obtain sharp en- 
closures for the simple zeros of nonlinear equations. First, for theoretical 
considerations, we present the following lemmas. 


Lemma 1. [17, 1] For real numbers a and 0 and interval numbers x and y, 
we have 

(i) w(ax + by) = |alw(x) + |b|w(y), 

(ii) w(xy) < |x|w(y) + lylw(x). 


Lemma 2. [17] Every nested sequence {x;} converges and has the limit 
Np Xk: 


Lemma 3. [17] If F is a natural interval extension of the real-valued rational 
function f with F(x) defined for x C xo, where x and xp are intervals, then 
there exists a constant DL such that 


w(F(x)) < Lw(x). 


Now we introduce the interval extension of (2) as follows: 
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Yn = N(x,) a) Xn; 
Yn = R(Xn, Yn) NXn, (4) 
Xn+1 = S(Xn, Yn, Zn) a Xn, 


where 


in which H is the interval extension of the continuous rational function h. 


Now we are ready to present the theoretical analysis of the proposed 
method (4). 


Theorem 2. Assume that f : D C R > R is continuously differentiable and 
that 0 ¢ F’(xo) for a given xo C D. If xo contains a zero 2* of f(x), then so 
do all x, for k = 1,2,..., defined by (4). Furthermore, the intervals x, form 
a nested sequence converging to x*. 


Proof. Using the Taylor expansion around x € xo, we have 


O= f(a") = f(x) + (@* — 2) f(r), 


for some €; between x and a*. Because f’(€,) 4 0, we obtain 


~_,_ F(z) 
= f'(E1)’ 8) 
which f'(&1) E F’ (xo) yields 
eee £6 IO 


Since x € xo is arbitrary, so in particular for « = m(xo), and taking into 
account that x* € x9, we obtain 


a" €E {m(xo) _ a xo = N(xo) (\xo = Yo. 
Now again using the Taylor theorem, for y € yo, we can write 


f(y) = f(a") + (y— 2*) f'(&), 


for some &) between y and a*. Since f’(€2) 4 0 and taking into account that 
f(a*) = 0, we get 
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. fy) 
ge=y- : (9) 
ata) 
As previously mentioned, method (2) is based on the King family. King [14] 
proposed the following formula for approximating f’ (yn): 


f(@n) + YF Yn) 
f(a) + BF Gn) 
with y = 6 — 2 to achieve a fourth-order of convergence. Let € and €2 be 


sufficiently close to # and y, respectively. Whereas in method (2), 6 = -4, 
and using (10), we have 


f'n) © fan) (10) 


2f(x) — 5f(y) 


Plea) Feo) =i 


; (11) 
Substituting (11) into (9) yields 


._,_ £) 2f(a) — fw) fo) i 


ac =yY => 


fié) ” 2F(@)—5f(y) FG) 


Indeed f’ (1) € F’(xo) and (12) holds for any x € xo and y € yo, in particular 
for « = m(xo9) and y = m(yo). So, we obtain 


2f(m(xo)) — f(m(yo)) f(m(yo)) = = 
2° { myo) — 3 ree) af ae)) Fag) f= RF) = 
Now for z € Zo, by the Taylor theorem, we have 
f(z) = f(a*) + (e — 2*) f' (Gs), (13) 
for some €3 between z and «*. Using the Taylor expansion for h(j) around 


zero with p = Het we get 


h(j) © h(O) + uh! (0). 


Since h(0) = 1 and h’(0) = 2, we obtain 


and so 

f(z) 

f(a)” 

Because z is arbitrary, we can assume that z and 2* are sufficiently close 
together and so f(z)h(u) = f(z). Now using (13), we obtain 


f(z)h(u) = f(z) +2 
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oe ty £2 
feel ee 4) 


Indeed f’(€3) € F’(z9) and (14) holds for any x € xo and z € Zo, in particular 
for z = m(xo) and z = m(zo). Therefore, since x* € x9, we obtain 


f(m(Zo)) 
F’ (zo) 


ne {m(ao) H (io) 


\ M xo = S(Xo, Yo, Zo) M Xo = Xj. 


By continuing this process, we see that 
x* E Xz, k=0,1,.... (15) 


Now by formula (4), it is obvious that x,41 C xx for k = 0,1,..., which means 
that {x;,} is a nested sequence. By Lemma 2, this sequence is convergent to 
a=MP 1X. Since x* € x, for all k, then x* € a. On the other hand, m(z,,) 
is not contained in S(xn,Yn,Zn) unless f(m(z,)) = 0. Since m(z,) € Zn C Xn, 
we conclude that w(x,+41) < w(x,). Therefore a = 2*. 


Note that procedure (4) stops when some stopping criteria are fulfilled, 
such as w(x,) < € for a tolerance € or Xp 41 = Xp. The computational scheme 
of the proposed interval method (4) for enclosing the simple roots of a given 
nonlinear equation f(z) = 0 is presented in Algorithm 1. 


Algorithm 1 The new interval method (4) for enclosing roots of nonlinear 
equation f(x) = 0 


1: procedure INTERVAL ROOT-FINDING(f, Xo, tol) 
2: n= 0; 

3 while w(x,,) >= tol do 

4 Compute N(x,,) from (5); 

5 Yn =intersect(N(x,),Xn); 

6: Compute R(x, yn) from (6); 

7: Z, =intersect(R(Xn, Yn), Xn); 

8 Compute S(Xn,¥n,Zn) from (7); 

9: Xn41 =intersect(S(Xn,Yn,Zn),Xn)} 

10: n=n+1; 


11: end while 
12: return x, 


13: end procedure 


Theorem 3. Suppose that f : DC R - R is continuously differentiable and 
that 0 ¢ F’(x9) for a given x9 C D. If a* € xo, then x; contains a unique 
root of f(x), for k = 0,1,.... Furthermore, if S(xz, yz,2%) A x~ = 0 for some 
k, then f(x) £0 for all x € xo. 
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Proof. Let x* € x9. By Theorem 2, we conclude that «* € x, for all k, which 
is unique because 0 ¢ F’(x;,) C F’ (xo). 

Now suppose S(xx,yxz,Z2%) 1 x~ = @ for some k, but x* € xo is a root of 
f(x), so by Theorem 2, we conclude that 2* € x, for all n. Particularly 
xv* © Xp41 = §(Xk, Ye, Zk) Xx, which is a contradiction. 


Theorem 3 is in the category of verification methods. By verifying its 
assumptions with the aid of a computer, we can detect when a certain interval 
does not contain a root. 


Theorem 4. Let f : DC R > R be continuously differentiable and have 
a simple zero x* in xo. If 0 ¢ F’(xo), then the interval method (4) has an 
eight-order of convergence, that is, there exists a positive real number C’ such 
that 


w(Xn41) < C(w(xn))*. 


Proof. Since Xn41 © S(Xn,Yn,Zn) SO W(Xn41) < W(S(Xn,Yn,Zn)). By the 
mean value theorem, we can write 


f(m(2n)) = f(m)(m@n) — @*), 
for some 7 between m(z,,) and x*. Thus we get 


f'(m) (a(n) = 2°) 


S(Xn,Yn; Zn) = m(Zn) _ H(fin) F’ (Zn) 


(16) 
Therefore, from (16) and Lemma 1, we obtain 


W(S(Kn,¥nsZn)) < |H(fin)||m(Zn) — IF mw ees) + 
W(H(fin))|m(2n) — 2*||F"(m)|le7egl- (17) 
Because x* € Zn, it is obvious that 
|m(Z,) — x*| < w(Zn). (18) 
On the other hand, we can write 


2f(m(Xn)) = f(m(yn)) f(mn)) 
2f(m(Xn)) — 5f(m(yn)) Fn) - 


Using the mean value theorem, we have 


f(m(Yn)) = f"(m2)(m(yn)-2"), and sf (m(Xn)) = f"(m3)(m(Kn) — 2") 


Zn © m(Yn) (19) 


for some 72 between m(y,) and xz* and some 73 between m(x,,) and 2*. 
Because 


Im(¥n) —2*| <wl¥n) < w(K) and — |m(X,) —2*|<wkn), (21) 
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so using (20) and (21), we obtain 
|f(m(¥n))| < LF (m2) Wn), (22) 


2f(m%n)) — f(m(Yn))| < 21 F(mEn))| + |F(m(yn))| (23) 
= 2| f"(ns)||m(%n) — «*| + [f"(n2)|lm(yn) — 2" | 
< | f"(ns)\w(Xn) + Lf" (n2)/w(%n) < Crw&n), 


where C, is an upper bound for 2] f’(m3)| + |f’(m2)|.. On the other hand, 
Theorem 1 yields 


w(Yn) < Co(w(%n))?, (24) 
for a positive constant C2. So by (22) and (24), we obtain 
|f(m(¥n))| < C3(w(%n))?, (25) 


in which C3 is an upper bound for C2] f’(72)|. Using Lemma 3, we have 


) < Cyw(x,). (26) 


Now from (19) and Lemma 1, we can write 


l2f(m(xn)) — Flon(¥n))| 1 
S [ap m(x,)) Sf tm(yn) yo mm) re 


Moreover, using (23), (25), and (26) yields 


w(Zn) ); 


W(an) < C5 (w(Xn))*, (27) 


where Cs is an upper bound for BAGG NS EGae By Lemma 3 and (27), 


we get 
1 


Therefore, from (18), (27) and (28), we have 
~ I) pl 1 
|H (fin )||m(n) — 2" || f (mle y) < Ce(w(Xn))®, (29) 


in which Cg is an upper bound for C?|f’(7)||H(fn)|.. Now by Lemma 3, 
there exists a positive constant C7 such that 


Using Lemmas 1 and 3 and (27), we obtain 
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gp Fl@n)  _ WFn)) Calan) 

wllin) = WF)? [Flanlen))] = [FtanGen)) 
CsCg 


where Cx is a positive constant and C5 is an upper bound for TFmGeT From 
(30) and (31), we can write 


) = < Co(w(xn))*, (31) 


w(H(jin)) < Cro(w(xn))*, (32) 


in which Cio = C7Co. Using (18), (27), and (32), we obtain 


w(H(fin))|m(Zn) — ellie < Ci(w%n))*, (33) 


where C, is an upper bound for C5C jo] f’(71)|| Fas) |. Finally, since w(Xn4i1) < 
w(S(Xn,;Yn,Zn)), by (17), (29), and (33), we conclude that w(xn+i) < 
C(w(x,))®, where C = Cg + Cy. 


As one can see, the new interval method (4) with three-step has an eight- 
order of convergence, while some other interval methods with the same num- 
ber of steps have a lower order of convergence; for some of them, see [1, 21]. 


3 Test problems 


In this section, we give some numerical examples to illustrate the performance 
of the new approach proposed in Section 2. The new method is compared 
with the interval Newton method, interval Ostrowski method, and interval 
modified Ostrowski method. In all examples, the procedures are stopped 
when w(x) < 10~1®. We utilize INTLAB [22] to compute the verified results 
on the computer. We study the following examples: 


1 
fi(a) = aresin(x® —1)—5a+1, aj ~0.5948109683983692, 
fo(x) =In(a?+242)—a241, xh © 4.1525907367571583, 
fa(x) =a? —e 3442, wf  0.25753028543986079, 
fu(x) =arctan(z)+2—8, 2% © 6.58002470991429699, 
fs(a) =a -1/2, ee=1. 


The first two examples are taken from [2] and the latest is taken from [15]. 
For all examples, we use rational function h as follows: 


In Figure 1, one can see the graphs of five functions f1, f2, fs, fa, and fs, 
respectively, over the initial intervals x} = [0.4, 1], x2 = [3.5, 5], x3 = [0.1, 2], 
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fi (2) § aresin(w? —1)—3a2+1 fe(x) =In(a?+24+2)—a24+1 _fs(x) = a? —e* —344+2 


5 
0.4 
ie} 0 
0.2 
0.5 5 
ie} 
-0.2 -1 -10 
0.4 0.6 0.8 1 3.5 4 4.5 5 0.5 1 1.5 2 
; = aretante _8 fs(x) =ax—1/a 
A fa(z) = arctan(x) + a a 
ie) 
2 
-0.5 
ie} 
-1 
2 =15 
5 6 rs 8 9 0.6 0.8 1 1.2 


Figure 1: Graphs of functions f1, fo, fs, fa, fs 


x} = [5,9], and x = [0.5,1.2]. By this figure, in addition to obtaining 
an intuitive view of the functions, we can understand the behavior of the 
functions for computing the following parameter: 


px = maxx |f()| (34) 


In the tables below, one can see the results obtained by implementing the 
interval Newton method, the interval Ostrowski method, the interval mod- 
ified Ostrowski method, and the new method (4) introduced in this paper. 
The third and fourth columns of the tables show, respectively, the tolerance 


parameters 6, = SS: and px introduced by (34). Note that in some 
tables, mark ”—” in the last step shows that the method fails in solving the 
problem. 


As the first example for f,(x) = arcsin(x? — 1) — $x +1, we present the 


obtained results by the mentioned methods in Tables 1-4. The presented 
results in these tables show that the new method (4) achieves the desired 
result with less number of iterations and higher accuracy. Also, the interval 
Ostrowski method fails in solving the problem. 

Tables 5-8 show the results obtained by executing different methods for 
enclosing the root of fo(x) = In(a?+2+2)—2+4+1. It can be seen that the 
new method (4) successes in the least number of iterations. Also, the interval 
Ostrowski method fails in solving the problem. 

For the third function f3(x) = 2? — e” — 3x + 2, the reported values in 
Tables 9-12 show that only the new method (4) successes in getting the result 
and the other methods fail. 


Iran. j. numer. anal. optim., Vol. 13, No. 1, pp 102-120 


Dehghani-Madiseh 


112 


Tables 13-16 display the results obtained by executing four methods for 


enclosing the root of f4(”) = arctan(x) + x — 8. As one can see, the interval 
Ostrowski method has failed to obtain a result, and the new approach gives 
better results than the other methods. 


The results of different methods for obtaining appropriate enclosures for 


the positive root of f5(x) = « — 1/x have been displayed in Tables 17—20. 
The interval Newton method does not yield any result. Whereas the new 
approach yields the exact root of the function. 


Table 1: Results of the interval Newton method for f1(a) = arcsin(x? — 1) — $a +1 


ko Xk OK Pk 
1 [0.40000000000000002, 0.66396313641487115] 2.64 x 10-7 747 x 10°? 
2 [0.56560254826011236, 0.66396313641487115] 9.84 x 107? 7.47 x 10~? 
3 [0.59018815218397114, 0.59856980551945871] 8.38 x 1073 3.98 x 1073 
4  [0.59480310218157917, 0.59481912020532601] 1.60 x 10-® 8.63 x 10~ 
5 [0.59481096839332148, 0.59481096840342751] 1.011071! 5.36 x 10- 
6 [0.59481096839836900, 0.59481096839836911] 2.22x10-1® 0 
7  {0.59481096839836911, 0.59481096839836911] 0 0 


Table 2: Results of the interval Ostrowski method for f(x) = arcsin(x? — 1) — $a +1 


ko xR Ok Pk 

1 [0.54158214865149934, 0.63394129754193074] 9.24 x 10-7 4.19 x 10-2 
2 [0.59477478728793232, 0.59485799844400755] 8.32 x 10-5 4.98 x 107° 
3 0.59481096839836756, 0.59481096839837055] 3.11 x 107! 1.33 x 107-15 
4 


Table 3: Results of the interval modified Ostrowski method for f1(2) = arcsin(«? — 1) — 


$a +1 
ko Xp Ok Pk 
1 [0.58885410911304559, 0.59936304066316770] 1.05 x 10-2, 4.83 x 10-3 
2 [0.59481096839549719, 0.59481096840132608] 5.83 x 10-% =3.13 x 10-7 
3 [0.59481096839836922, 0.59481096839836933] 2.22 x 10-16 =1.11 x 10716 
4  [0.59481096839836933, 0.59481096839836933] 0 1.11 x 10716 
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Table 4: Results of the new method (4) for fi(x) = arcsin(x? — 1) — $a +1 


Xk Ok Pk 


k 
1 
2 
3 


[0.58015286826057066, 0.60890961953980971] 2.88x 10-2 1.50 x 107? 
[0.59481096839720404, 0.59481096839958292] 2.38x 107 1.29x10- 
[0.59481096839836911, 0.59481096839836911] 0 0 


Table 5: Results of the interval Newton method for fo(a) =In(#? +a2+2)—-a2+1 


Xf Ok Pk 


oR WN EH! 


[4.09482718955130400, 4.17132082850488750] 1.83x 10-2 3.47 x 107? 
[4.1523 1696283340760, 4.15292802943720400] 1.47x 1074 2.03 x 107+ 
[4.15259073289156170, 4.15259074074274000] 1.90 10-9 2.40 x 107° 
[ 
[ 


4.15259073675715750, 4.15259073675715840] 4.28x 10716 4.44 x 10716 
4.15259073675715840, 4.15259073675715840] 0 0 


Table 6: Results of the interval Ostrowski method for fo(a) =In(a? +2+2)—-a+4+1 


Xk OK Pk 


WNrF| a 


[4.14427225093898070, 4.15515943057456380] 2.62 107% 5.01 x 107% 
[4.15259073560489430, 4.15259073791874480] 5.57x 10718 6.10 x 10710 
[4.15259073675715840, 4.15259073675715840] 0 0 


Table 7: Results of the interval modified Ostrowski method for f2(«) = n(x? + x + 2) — 


xa+1 
k Xp Ok Pk 
1 [4.15136705154255560, 4.15297239536206850] 3.87 x 10-4 7.37 x 1077 
2 (4. 15259073675715750, 4.15259073675715840] 4.28 x 10-16 4.44 x 10716 
3 — — 


Table 8: Results of the new method (4) for f(x) = In(a? +2+2)—a2+1 


ko Xp dk Pk 
1 [4.15167922809522590, 4.15321948581378480] 3.71 x 10-4 5.49 x 10-7 
[4.15259073675715840, 4.15259073675715840] 0 0 
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Table 9: Results of the interval Newton method for f3(x) = 2? — e” — 32 +2 
ko Xp Ok Pk 
1 — [0.10000000000000001, 0.76487534371627797] 6.65 x 1071 1.86 
2  [0.17953909948997981, 0.30082399330312792] 1.21x10-! 2.97x 107! 
3 [0.25663052647850410, 0.25844642836458781] 1.82 10-% 3.46 x 1073 
4 [0.25753027894621072, 0.25753029191301735] 1.30x 10-8 2.45 x 1078 
5 (0. 25753028543986067, 0.25753028543986073] 1.11 x 10716 4.44 x 10716 
6 — — 


Table 10: Results of the interval Ostrowski method for f3(x) = 2? — e” — 3a +2 


Xk OK Pk 


BWM FR] =a 


[0.10000000000000001, 0.31655239623745746] 2.17x10-' 6.05 x 1071 
[0.25752321108442017, 0.25753842849505237] 1.52 x10-° 3.08 x 10-° 
[0.25753028543986073, 0.25753028543986078] 1.11 x10-'S 0 


Table 11: Results of the interval modified Ostrowski method for f3(x) = x? —e® —3x42 


ko Xp Sk Pk 
1 [0.24154741311026207, 2.00000000000000000] 8.79 x 10~' 7.39 
2  [0.25749104640972659, 0.39675078835778121] 1.39 107' 5.20x 107! 
3 (0. 25753043640242368, 0.25753384076872499] 3.40x10-° 1.34x10-° 
4 = _ 

Table 12: Results of the new method (4) for f3(x) = x? — e? — 3x42 
ko Xp Ok Pk 
1 [0.22110828457567316, 0.27623770073133980] 5.51 x 107? 1.38 x 107! 
2  [0.25753028543982470, 0.25753028543989787] 7.321074 1.40x 107! 
3 [0.25753028543986078, 0.25753028543986078] 0 0 
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Table 13: Results of the interval Newton method for f4(x) = arctan(x) + « — 8 


Xk Sk Pk 
[6.5762681889199976482, 6.5869858860385530619] 1.62 x 10-3 7.11 x 10-3 
[6.5800246452848929479, 6.5800247578416417582] 1.71 x 10-8 6.60 x 1078 
[6.5800247099142961105, 6.5800247099142978868] 2.69 x 10-16 §=1.77 x 10715 
[6.5800247099142969986, 6.5800247099142969986] 0 0 


BRwWn eH) 2 


Table 14: Results of the interval Ostrowski method for f4(x) = arctan(x) + « — 8 


ko XxX, Ok Pk 
1 [6.5799958235806119689, 6.5800370828300822623] 6.27 x 10-6 2.95 x 107° 
2  [6.5800247099142961105, 6.5800247099142969986] 1.34 x 10-18 8.88 x 10716 
3 —= —, — 


Table 15: Results of the interval modified Ostrowski method for f4(x) = arctan(x)+a2—8 


k Xp Ok Pk 
1 [6.5800246462005800296, 6.5800248588084278012] 3.23 x 10-8 1.52 x 10-" 
2  [6.5800247099142969986, 6.5800247099142969986] 0 0 


Table 16: Results of the new method (4) for f4(x) = arctan(x) + #— 8 


ko xXx Ok Pk 
1 [6.5800247087713694683, 6.5800247 104028359857] 2.47x 10-19 1.16 x 1079 
2 [6.5800247099142969986, 6.5800247099142969986] 0 0 


Table 17: Results of the interval Newton method for f5(x) = x — 1/a 


k xp OK pk 
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Table 18: Results of the interval Ostrowski method for f5(x) = x — 1/a 


k Xk Ok Pk 

1 [0.99046958119024919309, 1.0128785276723828446] 2.21x 10-2 2.55 x 107? 
2 [0.9999999856310709534, 1.000000014532778092] 2.89x 10-8 2.90 x 10-8 
3 [1.000000000000000000, 1.000000000000000222] 2.22 x 10716 = 3.33 x 10716 
4 — [1.000000000000000222, 1.000000000000000222] 0 3.33 x 10716 


Table 19: Results of the interval modified Ostrowski method for fs5(x) = # — 1/x 


ko Xp 5k Pk 
1 [0.99900511706023975567, 1.0007695812111181421]  1.76x 10-7 1.99 x 10-3 
2 — [1.000000000000000000, 1.000000000000000222] 2.22 x 10718 3.33 x 10716 
3 [1.000000000000000222, 1.000000000000000222] 0 3.33 x 10716 
Table 20: Results of the new method (4) for fs5(x) = a — 1/a 
k Xk Op Pk 
1 [0.99968995513425429333, 1.0004281041560696419] 7.37x 10-4 8.56 x 10-4 
2 [1.0000000000000000000, 1.000000000000000222] 2.22 x 10718 3.33 x 10716 
3 1 0 0 


4 Enclosing the roots using extended interval arithmetic 


In this section, we want to introduce a technique that can find all the roots 
of a nonlinear equation f(2) = 0 located in a wide initial interval. Many 
root-finding methods in floating-point arithmetic can only find one root of 
the function in a given interval. Our technique is based on combining the 
new method (4) introduced in this paper and the extended interval arithmetic 
[12, 17]. 

For a continuously differentiable function f(x), if xo contains more than 
one zero of f(x), then o € F’(xo), and the discussed theorems in Section 2 
will not be applicable. Using the extended interval arithmetic, this problem 
can be handled. As said in [17], the definition of interval division can be 
extended as follows: 


[a, b]/le, d] = [a, ](1/Te, dl), 
where 

1/[c,d] = {1/y = y € [e, d]}. 
If 0 ¢ [c,d], the we are using the ordinary interval arithmetic. If 0 € [c,d], 
leaving aside the case c = d = O, then the extended interval arithmetic 
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- f(a) = sinh(x) — 2” tan(z) : 


L Ll L L 
aa 0.5 0 05 1 15 


Figure 2: Graph of function f(x) = sinh(x) — x? tan(« 


specifies the following cases: 


[1/d, +00) ifc=0<d, 
1/[e, d] = 4 (—00,1/c] U[1/d,too) ife<0<d, 
(—00, 1/¢ ife<d=0. 


Now if the initial interval xo is such that 0 € F’(xo), then the quotient 
tee) in (5) splits into two unbounded intervals. Thereafter intersecting 
N(xg) with the finite interval x9 yields two disjoint intervals y,, and yjpo. 
First, for yi1, if 0 ¢ F’(yi1), then we take y1, as the initial point for the new 
method (4), otherwise again by computing N(yi11) and then intersecting it 
with y11, we obtain two other intervals. By repeating this process, we find 
some intervals that contain a simple zero of f(a) and F’ over them does not 
contain zero. The process for yj is similar. Considering these intervals as 
initial points for the new method (4), we find all roots of f(x) on the initial 
interval xp. A similar idea previously has been used for the interval Newton 
method; see [17]. 

For an example, we consider f(a) = sinh(x) — x? tan(x) on the initial 
interval xo = [—1, 1.5]. The graph of this function on x9 = [—1, 1.5] is shown 
in Figure 2. We have 0 € F’(x) = 10?[—4.9097, 0.3056]. Using the extended 
interval arithmetic, we obtain 


f(mid(xo)) 
F’ (xo) 
= (—o0, 0.24225490053166] U [0.25048201511344, +00). 


N(xo) = m(xo) — 


Intersecting N(x) with xo, we get 


yi = [—1, 0.24225490053166] U [0.25048201511344, 1.5). 
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Indeed F’({—1,0.24225490053166]) and F’((0.25048201511344, 1.5]) contain 
zero, too. We repeat the above process by putting xg = [—1, 0.24225490053166] 
and xo = [0.25048201511344, 1.5], separately. Doing this work several times, 
we obtain three appropriate intervals, and then we apply the new method 
(4) on these intervals. The obtained results are shown in Table 21. 

As one can see, in a few iterations, all three roots of f(x) = sinh(x) — 
x* tan(x) in x9 = [—1,1.5] have been enclosed with sharp bounds and high 
accuracy. 


Table 21: Results of the new technique in Section 4 for f(x) = sinh(x) — x? tan(«) 


Xk OK Pk 

X11=[0.87539095698495750, 0.93264721412667118] 5.73 x 10-? 9.89 x 10-? 
X12=[0.90196399818943263, 0.90196401144268923] 1.33 x 1078 2.08 x 10-8 
X13=[0.90196400520858955, 0.90196400520858966] 2.22x 10-16 4.44 x 10-16 


X21=[-0.00000014204496225, 0.00000008340020063} 2.25 x 10-7 ~=1.42 x 10-7 
X2o=107°°[-0.20045735325692, 0.46773382426614] 6.68 x 10-5! 4.68 x 1075! 
X31 =(-0.90414914681585001, -0.90027520645356984] 3.87x 10-3 6.51 x 10-8 
X32=[-0.90196400520858988, -0.90196400520858899] 8.88 x 107'® 1.33 x 10748 


Nr Nr Wwne|] w= 


5 Concluding remarks 


In this work, a new family of numerical methods for enclosing the simple roots 
of the nonlinear equations was proposed. We showed that the new methods 
have an eight-order of convergence and also that the convergence analysis of 
the methods was studied. Some numerical examples were presented to show 
the feasibility and effectiveness of the new method proposed in Section 2. 
Also, we proposed a technique based on combining the new method (4) with 
the extended interval arithmetic to find all the roots of a nonlinear equation 
located in an initial interval. Finally, a numerical example for testing this 
technique was presented. 
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